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Abstract 

Marginal LTB models with corrections from loop quantum gravity have recently 
been studied with an emphasis on potential singularity resolution. This paper cor- 
roborates and extends the analysis in two regards: (i) the whole class of LTB models, 
including non-marginal ones, is considered, and (ii) an alternative procedure to de- 
rive anomaly-free models is presented which first implements anomaly-freedom in 
spherical symmetry and then the LTB conditions rather than the other way around. 
While the two methods give slightly different equations of motion, not altogether sur- 
prisingly given the ubiquitous sprawl of quantization ambiguities, final conclusions 
remain unchanged: Compared to quantizations of homogeneous models, bounces 
seem to appear less easily in inhomogeneous situations, and even the existence of ho- 
mogeneous solutions as special cases in inhomogeneous models may be precluded by 
quantum effects. However, compared to marginal models, bouncing solutions seem 
more likely with non-marginal models. 



1 Introduction 

Quantum gravity changes the structure and dynamics of space-time on small distance 
scales, which should have implications for the final stages of matter collapse. An interesting 
class of models to shed light on this issue is given by Lemaitre-Tolman-Bondi (LTB) space- 
times, which are inhomogeneous but do not show too much complexity. Classically, these 
models describe collapsing dust balls, containing Friedmann-Robertson-Walker solutions 
as special cases. They thus provide an interesting extension of models beyond homogeneity, 
an extension which is particularly important to understand in the case of quantum gravity. 
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Loop quantum gravity implies characteristic correction terms in the Hamiltonian con- 
straint of gravity and matter. Once a combination of corrections keeping the algebra of 
constraints anomaly-free has been found, a canonical analysis of quantum gravitational 
collapse becomes possible. Finding such equations is not easy, making the investigation of 
implications from quantum gravity corrections highly restricted. Quantum geometry cor- 
rections result from an underlying spatial discreteness, which may cast doubt on whether 
they can leave the theory covariant. There are arguments at the level of the full theory 
[lj stating that quantum operators may be anomaly-free, but it remains unknown how 
to descend from this statement to anomaly-free effective space-time geometries. A phe- 
nomenological approach has thus been followed to investigate possible geometrical and 
physical effects of diverse corrections. Here, one inserts expected corrections in the classi- 
cal constraints, suitably parameterized to reflect quantization ambiguities, and evaluates 
conditions under which the corrected constraints remain first class. As several articles 
have by now shown, it is indeed possible to have anomaly-freedom even in the presence of 
quantum corrections resulting from spatial discreteness [2J[3JII]. 

Marginal models, which are a subclass of general LTB models, have been analyzed in 
this spirit in Ref. |3j. This has resulted in consistent deformations which implement some 
types of quantum corrections without spoiling general covariance, and made possible an 
initial analysis of implications regarding effective pictures of collapse singularities. (At 
the fundamental level of dynamical difference equations in a loop quantization, spherically 
symmetric models are singularity- free [5] as are homogeneous models [6j[7J[E].) It turned 
out that there is no clear generic avoidance of either space-like or null singularities by an 
obvious mechanism, in contrast to several homogeneous models of loop quantum cosmol- 
ogy [9] where phenomenological mechanisms such as bounces could be found easily. While 
this outcome is not entirely unexpected given the types of corrections analyzed in the 
marginal case, it does show that further analysis is required. Marginal models, after all, 
provide spatially flat Friedmann-Robertson-Walker models in the homogeneous limiting 
case which give rise to phenomenological singularity avoidance in their loop quantization 
(including a positive matter potential) only with holonomy corrections [TQHH], which were 
not fully included in [3J due to technical complications. It is thus natural to extend the 
constructions to non-marginal models which would provide a homogeneous model of posi- 
tive (as well as negative) spatial curvature as limit. In that case, loop quantum cosmology 
can give rise to phenomenological singularity resolution even in the presence of inverse 
triad corrections alone [12], which in inhomogeneous situations are easier to control than 
holonomy corrections. If the behavior seen in homogeneous models should be generic and 
apply also to inhomogeneous situations, loop quantized non-marginal LTB models must 
give rise to singularity resolution more easily than marginal ones. 

To find anomaly-free versions of non-marginal LTB models including inverse triad cor- 
rections from loop quantum gravity, we will follow two derivations. First, we will extend 
the methods of [3J where constraints already incorporating the LTB reduction of metric 
components are made anomaly-free by consistency conditions between correction func- 
tions. Secondly, we will derive a general anomaly-free system of spherically symmetric 
constraints, on which we then apply the LTB reduction in a second step. As we will show, 
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the two steps of LTB reduction and deriving consistency conditions almost commute: in 
the end, we obtain consistent equations of motion of similar structure, although they do 
differ by some terms. This outcome considerably supports the constructions of [3]. 

Using these consistent equations, gravitational collapse can be analyzed. We are specif- 
ically interested here in the possibility of a turn-around of the collapse, or a bounce, in the 
corrected equations, which are suggested to exist by models where homogeneous interiors 
have been matched, Oppenheimer-Snyder-style, to spherically symmetric exteriors [T3] . 
Also here, as in the marginal case but in contrast to homogeneous models, we do not find 
a clear indication for singularity resolution, although several extra terms do seem to make 
a bounce more likely. As in the marginal case, this part of the result is not conclusive 
since not all corrections have been included and no complete analysis has been performed. 
Our results thus do not mean that there is no bounce in these inhomogeneous models. 
But they do show that an outright treatment of inhomogeneous models is different from 
matching homogeneous results. In fact, we also confirm the observation of [3 J that quan- 
tum corrections of the type studied here prevent the existence of an exact homogeneous 
limit. "Effective" homogeneous geometries thus have to be taken with care, but consistent 
relationships with inhomogeneous ones do provide insights in their structure. 



2 Classical equations 

Non-marginal LTB models O [T3J, |TB] have a space-time metric given by 

ds 2 = -dt 2 + —, dx 2 + R 2 dn 2 (i) 

1 + k(x) 

with dQ 2 = dfi 2 + sin 2 ddtp 2 and where k ^ is a function of the radial coordinate x. (The 
limiting case k = is that of marginal models.) The function R(t,x) can depend on both 
time and the radial coordinate, but not on the angular coordinates to leave the metric 
spherically symmetric. It is easy to see that positively curved Friedmann-Robert son- 
Walker models with scale factor a(t) are obtained for k(x) = — x 2 and R(t,x) = a(t)x. 

For an application of loop quantization we use densitized triads instead of the spatial 
metric components, whose conjugate momenta are given in terms of the Ashtekar-Barbero 
connection and extrinsic curvature components. Written as a densitized vector field taking 
values in su(2) with basis r», the spherically symmetric densitized triad is 

E = E x (x)r 3 sintf— + {E\x)r l + E 2 {x)r 2 ) sintf— + (E\x)t 2 - E 2 {x)r l )— . 

ox ov dip 

Similarly the Ashtekar connection A l a = T l a + 7-fQ, where and K l a are the components 
of spin connection and extrinsic curvature respectively and 7 is the Barbero-Immirzi pa- 
rameter [T71 IE] , reads 

A = A x (x)r 3 dx + (Ai(x)ri + A 2 (x)T 2 )d-d 

+(Ai(x)t 2 — A 2 (x)ti) sin-iM^ + r 3 cos^dv? 
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Introducing the U(l)-gauge invariant quantities (E^) 2 = (E 1 ) 2 + (E 2 ) 2 and A 2 = A\ + A\ 
(see [19l [20j EI] for details) we have the symplectic structure 

{A x (x),E*(y)} = { 7 ^(x),2^(y)} = { V (x),P*(y)} = 2G 1 6(x,y) 

or more explicitly the Poisson bracket of functions / and g is 



{f,g} = 2G J da; (7 



Sf 5g 1 8f 5g Sf 5g 



SA X SE X 2 5K v 8Ev ! 5r]SP^ 
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Compared to metric variables, we have an extra field rj(x) with momentum 

= sin a = 4tr {{E 1 t x + E 2 r 2 )(A 2 r 1 - A x r 2 )) 

(with a defined as the angle between the internal directions of A- and ^-components), 
which plays the role of a U(l)-gauge angle in the spherically symmetric theory. (This 
gauge angle also determines the x-component of the spin connection T x = —rj', and thus 
enters the Ashtekar connection by A x = —r)' + jK x with an extrinsic curvature component 
K x .) 

In this situation, we have three constraints: the Gauss constraint 
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W = 2^/ teKm' + P*) (2) 



the vector constraint 



D^[N X ] = ^JdxN x (2E"K' V - ^A X (E*)> + UfP^j 

^JdxN x (2E*K' V - K X (E X Y + \'{{E X )' + pn] 



2G 

and the Hamiltonian constraint 



(3) 



H gTav [N] = ~ J dxN\E*\-h(KZE* + 2K V K X E X + (1 - T^E^ + 2V' ip E x ) (4) 

with r v = —(E x )' j2E v the gauge invariant angular component of the spin connection. 
Solving the Gauss constraint removes the pair (77, P v ) and reduces the vector constraint to 
the diffeomorphism constraint. After this step we can work with the canonical pairs 

{K x (x), E x {y)} = {K^x), 2E*(y)} = 2GS(x, y) . 

The relation to the usual spherically symmetric geometrodynamical variables 

{R(x),P R (y)} = {L(x),P L (y)} = GS(x,y) (5) 
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as used for example in [22j [23] can be obtained directly by comparing the spatial metric 

dq 2 = L 2 dx 2 + R 2 dn 2 = ^tdx 2 + \E*\dtf (6) 

\E X \ 

in each set of variables and making use of the equations of motion: 

L = E*'\E X \-* , R=\E X \*, 
P L = -K v \E x \ l * , P R =-sK x \E x \*-K ip E v '\E x \-* (7) 

where s = sgn(E x ). (The sign factor corresponds to the two possible orientations of a 
triad; we will mostly use s = +1 below.) 

Specializing the general spherically symmetric metric 

ds 2 = -N(t, xfdt 2 + L 2 (t, x)(dx + N x (t, x)dt) 2 + R 2 (t, x)dtt 2 

to the LTB form ([1]) requires a vanishing shift function N x = and lapse iV = 1 for 
comoving coordinates of the dust, and on using the first equation in (JTj) gives the non- 
marginal LTB condition 

2y/l + k(x)E* = (E x )' (8) 
in terms of triads. From this we can derive the spin connection component 

r ^ = = -v / i+^y (9) 



and its derivative r' = —k'{x)/2^/1 + k(x), which appear in the Hamiltonian constraint. 
Since the spin connection, unlike in the marginal case, is not a constant —1, the Hamilto- 
nian constraint is different from the marginal case: 

H*Z S IN] = ~ J dxN^lET 1 / 2 (k%E* + IK^KJ? - k(x)E* - ^== 

(10) 

If we solve the diffeomorphism constraint identically, which requires 2E Lp K'^—K x (E x )' = 
0, the LTB condition for triad variables gives rise to a condition 



K' v =y/l + K(x)K x (11) 

for the extrinsic curvature components. For a consistent LTB formulation, the two LTB 
conditions must be preserved by evolution generated by the constraint. This is indeed the 
case as can be seen from deriving Poisson brackets between the Hamiltonian constraint 
and each of the LTB conditions. For the Poisson bracket of the two LTB conditions, after 
smearing them with fields fi(x) and u(x), we get 

J dxu(x) ( y/i + k{x)k x -k)>[ d yv(y) (2 V / i + k(x)^ - (E x y 

2G [ dzy/1 + k(z)(hv)' . (12) 



This in general is non-vanishing, unlike in the marginal case where k is zero. (Although we 
will not follow this route here, we note that this will have an impact on implementing the 
LTB conditions at the state level, as done in [3] for the marginal case. Another complication 
for such a construction is the explicit ^-dependence of the LTB conditions, which makes 
their integrated version used as conditions on holonomies more complicated.) 

Equations of motion in this canonical formulation are derived using E x = {E x , H^.^}, 
with a similar equation for E v . With these we can first eliminate K x and K v , and finally 
E v using the non- marginal LTB condition to obtain an equation entirely in terms of E x . 
After replacing E x by R 2 we obtain 



2Ga/1 + K{X) 

which has to be equated to the matter part of the Hamiltonian for dust given by H^ nst = 
— ^qF'/ a/1 + k(x). (A more general canonical derivation of the gravity-dust system will 
be given in Sec. 14.21 ) Thus 

2RRR + R 2 R' - k'R - kR' = F' (14) 



is the equation of motion, in agreement with the spatial derivative of RR? = k(x)R + F(x) 
which is the equation obtained by solving Einstein's equation for the non-marginal case. 



3 Inverse triad corrections from loop quantum gravity 

We will now repeat the canonical analysis using a Hamiltonian constraint containing cor- 
rection functions as they are suggested by constraint operators in loop quantum gravity. 
Consistency will then require conditions for the possible terms, which show how quantum 
corrections can be realized in an anomaly-free way. We discuss here only inverse triad 
corrections which are easier to implement, and which already provide insights into one of 
the main classes of quantum geometry corrections. 

Inverse triad corrections arise from every Hamiltonian operator quantized by loop tech- 
niques, where inverse components of the densitized triad appear. Such corrections are 
directly related to spatial discreteness of quantum geometry since densitized triads as ba- 
sic variables are quantized to flux operators with discrete spectra containing zero Since 
such operators do not have densely defined inverses, no direct inverse operator is available. 
Instead, well-defined quantizations exist based on techniques introduced in [H[25], implying 
corrections to the classical inverse. 

In several symmetric models, inverse triad operators and the corrections they imply 
can be computed explicitly [26J. As an example, spherically symmetric models used in [3] 
give rise to a correction function of the form 

.( A)=a #±Mi#EM^ (15) 

JZ-p 
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Figure 1: The correction functions a (A) (solid line) and /(A) (dashed line) where A is 
taken relative to A» := ^ / y/2£p. 



where A is the size of an elementary plaquette in a discrete state underlying an LTB 
geometry. For corrections in inverse powers of E x , which is proportional to the area of 
a spherical orbit, the relevant operators give rise to a dependence on plaquette sizes on 
orbits. For a nearly spherical distribution of N(E X ) such plaquettes making up the whole 
orbit, we thus have A = E x /N. Since we refer only to the orbit size, corrections thus 
naturally depend on E x only but not on E v . (That this is required will later be shown 
independently when we use anomaly freedom to rule out that a could depend on E 9 .) Such 
a correction function then multiplies any classical appearance of (E x )~ 1 in a Hamiltonian 
operator. In particular, classical divergences of inverse factors of E x are cut off as one can 
see from the plot in Fig. [TJ Correspondingly, the dynamics given by such a Hamiltonian 
will change from quantum corrections. Classically, i.e. for £p — > 0, we have a(E x ) = 1, 
and this limit is approached for large E x . Also this behavior of the correction function is 
illustrated in Fig. [Q 

For a large number Jsf of discrete blocks, the scale A is reduced compared to E x and 
corrections from a can be significant even for large E x . Since typically N, in relation to 
the underlying state, is not a constant but would depend on the size E x , different kinds 
of behaviors can arise. This phenomenon of lattice refinement [27J [28] is important to 
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capture the full dynamics of quantum gravity and its elementary degrees of freedom. It is 
also crucial for realizing the correct scaling behavior of correction terms under changes of 
coordinates [29]. In this paper, we will mainly be looking at general implications in local 
equations of motion, where the value or behavior of N is not important. More detailed 
investigations could at some point provide restrictions on the possible form of J4(E X ), and 
thus give insights in the required behavior of discrete quantum gravity states. 



3.1 First version 

We turn to the case of inverse triad corrections called 'first version' in |3j , where only those 
terms in the Hamiltonian with explicit 1 / v/ji£*] dependence are corrected by a factor 
a(E x ). Starting as in the marginal case, we first assume that the classical expression for 
the spin connection can be used, but show that this is inconsistent. The Hamiltonian, now 
assuming E x > 0, is 



,Mx)E' p k'{x)\/E x ~ 



4JM = -l/d,*(o(^ + affI1 .VF-.(ff, vei vi + k(x) 

(16) 

As in the marginal case, it turns out that the correction in the Hamiltonian can lead to 
consistent LTB-type solutions only if we also change the LTB conditions by a correction 
function f{E x ): 

(E x )' = 2y/l + K (x)f(E x )E* and K' v = y/l + n{x)f{E x )K x . (17) 

These relations still solve the classical diffeomorphism constraint identically, which does 
not receive corrections in loop quantum gravity. The main consistency condition then is 
that Poisson brackets of the LTB conditions with the Hamiltonian constraint vanish. 

Each of the Poisson brackets gives a differential equation for f(E x ): For the LTB 
condition corresponding to the triad variables we obtain 

2E X ^L = f(l - a) . (18) 
dE x v ; ( J 

This is the same equation as found for the marginal case in [3J. For N = 1, for instance, 
the solution is given by 

Ci \fE^e~ a l 2 

fm = r ; y/vV ; (19) 

(yEF + y/E*-7e*/2) [VE^+ y /E x + 1 £ 2 p /2j 

for E x > 7 4/2 and 

cav^exp (-\a + \ arctan ( JWJ(y^/2 - WjX) 

f{E x ) = ^ V U- (20) 

^ + y/EF + 74/2 
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for E x < 7 4/2 (dashed curve in Fig. QJ. Here c x = 2^ and c 2 = 2 5 / 4 e 1 / 2 -^ 4 7~ 1/4 V /2 
are constants of integration fixed, respectively, by the condition lim^^oo f(E x ) — > 1 and 
by requiring that f(E x ) be continuous at E x = , ~f£p/2. 

It is clear that we obtain the same equation because there is no influence of k in the 
evaluation of the Poisson bracket, k not affecting the terms containing K x and K v in the 
Hamiltonian. However, the differential equation for f(E x ) obtained by demanding that 
the corrected LTB condition for extrinsic curvature components is also preserved in time 
gives the equation 



df VI + K,K,„K x f ay/T+KK ip K x f n'f n'a 



-2Vl + KK v K x VE x -4 r + - — * XJ - _ y XJ J= + ^^ = o (21) 

v dE x y/E^ \[E> 2\fE x 2\fE x 

which, due to the «/-terms, is different from and in fact inconsistent with equation (|18p 
obtained from the LTB condition for triads. 

3.2 Second version 

We now repeat the above procedure for the second version of the inverse triad corrections 
for which the Hamiltonian is corrected by a(E x ) in all terms and reads 

HiLiN] = ~ I dx N^2- (k%E« + 2K,K X E X - k(x)E> - ) • (22) 

^ J V E x V a/1 + k{x) J 

Again, the LTB conditions in the form 

(E x )' = 2yfl + K(x)g(E x )E* and K' v = y/l + k(x) 9 (E x )K x (23) 

solve the diffeomorphism constraint identically. However, for the same reasons as noted 
above the conditions cannot be consistent: since the terms containing K x and in the 
Hamiltonian do not involve k, the Poisson bracket for the triads gives the same result as 
in the marginal case where the differential equation was 

dg da 

a dE x=9 dE x ' (24) 
with the solution g(E x ) = a(E x ). (This solution is unique with the boundary condition 
imposing that g = 1 for large arguments.) 

For the Poisson bracket involving the condition on extrinsic curvature, on the other 
hand, the terms in the Hamiltonian involving k are important and we get a different 
result: 

-2 a/1 + K{x)aK v K x ^-^- + 2 a/1 + ^x)K v K x s/W^g '' <UJ 



dE x v v ' v dE x " 2^E X 

i n=r da k'ol „ , nr . s 
-K , VE x ——g + — = = (25) 



dE x " 2^E X 

As in the previous case the presence of k' terms spoils the consistency: Using a = g the 
first two terms cancel while the rest would require da/dE x = (a — 1)/2E X with a solution 
a = 1 + cVE x violating the classical limit at large arguments. 
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3.3 Inclusion of corrections in the spin connection 

A direct extension of the results from marginal to non-marginal models is thus impossible. 
Here we have an example for the information gained by a phenomeno logical treatment: 
LTB-type solutions require additional corrections to compensate inconsistencies seen so 
far. Such corrections may be more difficult to derive from a full Hamiltonian, but they 
follow directly from a phenomeno logical treatment. A successful consistent implementation 
thus provides feedback on the full theory: additional corrections required for consistency 
must eventually follow from the full theory just like the primary correction a followed from 
inverse triad operators. 

In particular, to resolve inconsistencies, we thus have to include further corrections 
in terms not affected yet, the chief candidate being the spin connection terms in the 
Hamiltonian constraint. They vanish in the marginal case, such that results from there 
do not provide much directions for more general models. Moreover, such terms are in fact 
more difficult to derive from a full Hamiltonian so that not much is known about their 
form. We will now look for corrections in the spin connection terms which are such that 
they combine with those already used to provide a consistent formulation. 

3.3.1 Implementation 

Classically we started with the expression = —{E x )'/2E Lp = — yT+ k(x). However, 
with quantum corrections to the Hamiltonian the LTB conditions are also corrected. For 
example for the second version the modified LTB condition (E x )' = 2 a/1 + ng{E x )E ip , 
implies that now the spin connection is 



We include the additional factor of g{E x ) in the Hamiltonian by replacing any occurrence 
of the classical spin connection r^ lass from (Q with g(E x )T c ^ ass . 
Furthermore, the derivative of the spin connection then is 



This introduces an explicit (E x )' in the LTB-reduced Hamiltonian, which would imply 
{H[N], H[M}} ^ even though the diffeomorphism constraint has been solved identically. 
The system would thus be anomalous. We are finally led to incorporate another correc- 
tion function multiplying T'^ in a function h(E x ) so that in the Hamiltonian r' is to be 
replaced with h(E x )(r^ ass y . The form of h(E x ) will be determined by the requirement of 
consistency. With all the possible corrections, the new Hamiltonian in the second version 
is 




(26) 




(27) 
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We now demand that this Hamiltonian Poisson commutes with the LTB conditions in 
( 123]) . Here we note that the terms containing the spin connection (and its derivative) in the 
Hamiltonian do not contain K x or K v and therefore in evaluations of the Poisson bracket 
with the first LTB condition there will be no changes. This leads to the same differential 
equation ( 1241) for g(E x ) as obtained earlier, implying g(E x ) = a(E x ). Evaluating the 
Poisson bracket of the Hamiltonian with the second condition and using the solution for 
g(E x ) gives a differential equation for h(E x ): 

&{ahy/W) _ a 2 



dE x 2\/~E x 

Here, we can thus have a consistent formulation for the non-marginal case for a suitable h 
by correcting the spin connection terms. 

We proceed in a similar manner for the first version of the inverse triad corrections. 
From the LTB condition for triads in (fTTj) we find that the spin connection would be 
replaced with f(E x )r^ ass , and the derivative of the spin connection receives a correction 
function l(E x ) in the form Z(£' a; )(r^ lass ) / . With these changes the Hamiltonian is 

1 r , _/aK 2 E v „ „ ,_ a E^ 



1 r /aK z E v , 

H^JN] = -—JdxN(—^ + 2K tp K x V& + 

(30) 



apEV KOipE* k'/v 7 ^ 



y/E* VE^ VT+k 

Evaluating the Poisson bracket of this with the LTB conditions and equating them to zero 
implies that the equation for f(E x ) is unchanged compared to (ITB1) . The other Poisson 
bracket then gives a differential equation for l(E x ): 



dE x 2\JE X 

Differential equations for the correction functions h and I are difficult to solve in general 
for given a and /. For a near center analysis done later we will need the lowest term in h 
and I in a power series expansion in x. Integrating ff29l) and fl3Tl) keeping only the lowest 
order term in a and / we find the solution 



p 



i 



M^0)=(-^) 2 ^ , l(x~0) = ^ w R* (32) 



valid near the center. 



3.3.2 Ambiguities 

The independence of a follows from the consideration of inverse triad operators in the 
full quantum theory. Although it is not determined uniquely in this way (see [301 [31] for a 
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discussion), the general shape of this correction function is known well. No such arguments 
exist for some other correction functions such as h, whose form is thus less clear. More 
ambiguities are thus expected to arise for it. 

A basic condition on multiplicative corrections is that they be scalar to preserve the 
transformation properties of corrected expressions under changing coordinates. Among the 
basic triad variables, E x is the only one free of a density weight and thus can appear in 
correction functions in an unrestricted way. The other component E v , on the other hand, 
is a density of weight one and would have to appear in combination with other densities 
to result in a scalar. If only triad components are considered for the dependence, the only 
other density would be (E x )'. Scalars made from these densities, such as E V /(E X )' are 
however unsuitable for corrections since they are not always finite. 

In the present situation, we use the function k for a non-marginal LTB model, which 
means that we have another density, k', at our disposal. Scalars of the form (E x )'/k' or 
E v /k' are well-defined for most functions k of interest, and can thus arise in corrections. 
This enlargement of the space of acceptable variables means that additional ambiguities 
can arise. In the next section we will see how several of these ambiguities can be fixed 
by an analysis of the constraint algebra. The equations of motion will remain structurally 
similar, so that we proceed for now with an analysis of the equations resulting from the 
treatment done so far. 

3.4 Equations of motion 

Given the consistency conditions between correction functions we can derive consistent 
equations of motion even without having explicit solutions for the differential equations 
( 1291) and (I3T1) . Once consistent constraints are available, the derivation follows the classical 
lines which we briefly illustrate first: The first order equation (in time) has already been 
worked out in flHj) , so that we can go on to the evolution equation. From E x = {E x , H^!^} 



we have 



E x 



(33) 




Similarly, using K v = {K^, H' 



} we obtain 




(34) 



Eliminating from the above two equations we obtain 




(35) 



which, using E x = R 2 , can be written as 



2RR + R 2 = k. 



(36) 
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Eqs. (JHj) and (136]) are automatically consistent, which can be seen explicitly by subtracting 
a time derivative of (TH|) from a space derivative of (1361) . 

The same procedure is then applied to constrained systems including consistent correc- 
tion terms. To get the first order equation in version two we use (1251) in the equations of 
motion E a = {E a , H^ ay } to solve for the extrinsic curvature components 

E x K E v 

K v = -^= and K x = —7= — —77— • (37) 
v 2aV& a^W E* v ; 

Using these along with the LTB condition E v = (E x )'/2\/l + ng in (1251) we rewrite the 
Hamiltonian in terms of E x only: 

jj 1 / (E x f(E x )' E X (E X )' (E X ) 2 (E X )' da 

#ot*v = — ^ — r~/= o,„ + 



grav 2G\ 8VTT^a 2 (E x )y 2 2 v / TT^a 2 v / ^ 2 v / TT^a 3 ^ d£* 

+ 1 - a 2 - to 2 - v ; - — I , 38 

where we have already used the condition g = a. When equated to the dust Hamiltonian 
after using E x = R 2 a first order equation in time ensues: 

R 2 R' 2RRR' 2RR 2 R' da 9 9x _. 

— — + - — —— + (1 — a — na )R — k ahR = F . 39 

or or a 6 dR 

To obtain the evolution equation we use 

JT 1 a — aq 2 — nag 2 + aK 2 
= {K„ H»} = ?L_» <L (40) 



2 

together with K v = E x j2a\[E x from (|37p . such that 

4£> a d£> v y y ; v 1 

With g = a and = i? 2 this becomes 

2 RR + f 1 _ 2^^) i? 2 = -a 2 (l - a 2 - K « 2 ) . (42) 



dlogfl, 

It is easy to see that this equation has the correct classical limit and (using (129]) ) is 
consistent with the first order equation. 

Proceeding in a similar manner for the first version of the inverse triad correction we 
find that the first order equation is 

aR 2 R' + 2RR'R + a(l - f - nf 2 )R' - n'flR = fF' (43) 

and the evolution equation 

2RR + aR 2 = -a(l - f - nf 2 ) . (44) 

Using (13"T1) along with (|24p one can verify explicitly that the first order and the second 
order equations are consistent with each other. 
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3.5 Effective density 

To interpret effects from correction terms it is often useful to formulate them in terms of 
effective densities rather than new terms in equations of motion. As in the marginal case, 
we use the Misner-Sharp mass defined by 

m =-(l- V A RV A R) (45) 

where A = (1, 2) corresponds to the t — r manifold. Writing the metric for spherical dust 
collapse as 

ds 2 = -At 2 + L 2 {t, x)dx 2 + R 2 {t } x)dn 2 (46) 
the equation for the Misner-Sharp mass becomes 

m =*(l + #-*!). (47) 

Classically R' 2 /L 2 — 1 + k and therefore the Misner-Sharp mass is m = (RR 2 — kR)/2. 
With an effective density defined in terms of the Misner-Sharp mass by 

TYl 

£eff = ±kGR*R> (48) 
we find that for the classical collapse it is 

pi 

C ^ SS = 8ttGR 2 R' ' (49) 

Here we have made use of the equation of motion R 2 R = kR + F which is obtained from 
the Hamiltonian constraint. This is in agreement with the 00 component of Einstein's 
equation, G 00 = SnGe{t,x) where e(t,x) is the dust density, implying e = F' /SiiGR 2 R' . 
Thus, classically the effective density, defined in terms of the Misner-Sharp mass, is the 
same as the dust density. Moreover, as expected, the expressions for the two are unchanged 
compared to those for the marginal case. 

We now proceed in the same way to find the effective density for the first version of 
the inverse triad correction. With the new LTB condition E 9 = (E x )' /2\/l + nf(E x ), the 
metric coefficient L = E^/^/E^ implies L = R'/\/l + nf(R) using E x = R 2 . Therefore 
the Misner-Sharp mass as defined in fj4T|) is now 

m ' = + # _(! + «)/»). (50) 



The corresponding effective density as implied by (1481) is 



J 1 ffF' / ,2 ,2 -2 x K'flR K'f 2 R\ . 

= SGff ( V + {a ~ l * 3f + 3K/ - R - V + -IT - -W-) (51) 
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where we have made use of (fTBj) after substituting for E x in terms of R, and of (jl3"j) . We 
note that this equation has the correct classical limit. 

For the second version of the inverse triad correction we have E v = (E x )' /2y/l + ng(E x ), 
implying L = R' + Kg(R) where as seen earlier g(R) = a(R). With this the Misner- 
Sharp mass is 

m // = |(l + J R 2 -(l + ^ 2 ). (52) 

Using the relation a' = R'da/dR (where the prime denotes derivative with respect to x) 
along with (1311 we find that the effective density is 

j j 1 ( a 2 F' 1 — a 2 , 9 9n 2 ,- 9 9 9 . da n'a 3 h K'a 2 \ 

(53) 

As in the marginal case, these effective densities imply that the near center expansion 
for the mass function F can have different behavior compared to the classical case, as 
discussed in Sec. [6j The matter contribution to the effective density, as given by the first 
terms of ( ISTT) and ( |53l) . is the same as in the marginal case. 

3.6 Quantum correction to the energy function k? 

Physically one would expect that the energy function n, which is related to the velocity of 
the dust cloud, should also receive corrections after including quantum effects. To derive 
those, we have to find an independent definition of k referring only to the constraints or 
evolution equations derived from them. One possibility, in the classical case, is to use ( |36|) 
whose right hand side only contains the energy function. Once brought into an analogous 
form, a corrected evolution equation can directly be used to read off a corrected energy 
function. Specifically for version one, where the evolution equation is given by (jS]), the 
effective energy function is 

Ki s = af 2 K-a(l-f 2 ) (54) 

while for version two, where the evolution equation is given by fj4*2l . the effective energy 
function becomes 

K ^ = a 4 K - a 2 (l - a 2 ) . (55) 

This correction in effect would imply that the near center expansion for k can be different 
for the quantum corrected equations as we will see when we come to the near center analysis 
below. 

4 Spherically symmetric constraints 

We have now several versions of consistent sets of equations of motion for non-marginal LTB 
models including inverse triad corrections as expected from loop quantum gravity. To make 
these equations consistent, we had to introduce several correction functions in different 
terms of the Hamiltonian constraint, which were then related to each other by consistency 
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conditions following from the requirement that the LTB conditions be preserved. Since 
there is some freedom in choosing the places and forms of corrections in the constraint as 
well as the LTB conditions, one may question how reliable such an analysis is regarding 
the structure of resulting equations of motion or implications for gravitational collapse. 

Before analyzing corrected equations of motion further, we now present an independent 
derivation which starts with a consistent set of corrected spherically symmetric constraints, 
and then implements the LTB reduction. As we will see, the structure of the resulting 
equations is nearly unchanged, while much less assumptions about different corrections are 
required. With these two procedures we thus demonstrate the robustness of consistently 
including corrections at a phenomenological level. Note that this would not have been 
possible had we chosen to fix the gauge generated by the Hamiltonian constraint in any 
way instead of dealing with the anomaly-issue head-on. 



4.1 Gravitational variables and constraints 

Quantum corrections due to inverse powers of the densitized triad are introduced in the 
Hamiltonian constraint (j4j) by functions which we initially assume to be of the general form 
a(E x , E v ) and a(E x , E 9 ) entering the Hamiltonian constraint as 

H^[N] = ~ J dxN(a\E x \- l 2K 2 ^+2saK v K x \E x \^ +a\E x \^E^ 

-a T \E x \-^lE tp + 2sa v T' tp \E x fi). (56) 

To account for possible corrections from the quantization of the spin connection, as sug- 
gested by the previous analysis, we have also introduced functions ar(E x , E v ) and ar(E x , E v ) 
in those terms. The only restriction so far is that we have the same a in the first and third 
term of the Hamiltonian constraint due to their common origin from the inverse \E X \~ 1 / 2 . 
The two main cases of interest here are a = 1 or a = a, corresponding to two versions of 
inverse triad corrections. 

We now proceed to make the corrected constraints anomaly-free before implementing 
LTB conditions. (For a similar analysis for dilaton gravity, see jl].) To ensure anomaly- 
freedom, we must determine conditions under which the system of constraints, including 
its corrections in the Hamiltonian constraint, remains first class. Computing the Poisson 
bracket { tfg av [M] , Fg av [N] } gives 

{ g av [M] , H^ v [N] } = D grav [aa r \ E x \{E*y 2 {MN' -NM')] 

- G giav [aa T \E x \(E^~ 2 (NM' - MN')rj'\ 

+ -4 / dx ( MN ' ~ NM'){aa r - aa r ) ^^- 
2G J Ef 

+ ^Jdx(MN'- NM') (a'ar - aaJ0 • (57) 

For a first class algebra the last two terms, which are not related to constraints, must 
vanish, providing conditions on the correction functions. The vanishing of the last term 
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implies ar oc a, upon which the third term gives «r oc a. Therefore to recover the classical 
limit we must have: 

cxr = ct , ar = ol . (58) 

Thus, anomaly freedom requires corrections to the spin connection terms to be only due to 
the inverse power of the densitized triad factors they contain. This may look contradictory 
to what we derived earlier, where additional correction functions such as h were needed. 
However, the previous case (where LTB conditions were used instead of the diffeomorphism 
constraint) implicitly makes h dependent on (E x )' as well: Comparing the correction terms 
we have 

and we can write, using g — a: 

k' dE x v 1 

Thus, to match the current equations the correction function h used earlier must depend on 
(E x )', which has a density weight. (Similar considerations apply to the correction function 
I.) As the expression demonstrates, this is made possible since in our earlier procedure we 
had the function k at our disposal in addition to the triad components. Its derivative n' 
provides an extra density, which can be combined with (E x )' to provide a scalar correction 
function. In the current setting, by contrast, we have not yet introduced any such function 
by LTB conditions, and so a possible dependence on (E x )' is more restricted. The new 
procedure of this section is clearly less ambiguous, while the final results will be very close. 
This again demonstrates the robustness. 

To continue with the analysis of anomaly-freedom, we compute Poisson brackets: 

__L J ^ N(N .y EV ^ lEriKiEV + 2s ^ KvK ^i + ^ Er i Ef 

-^ET L <E* + 2 S ^mi). (59) 

In the case a — a, 

{H^JNID^N*}} = -H^[N X N' - {d\oga/dE*)E*N{N*)'\ . 

The corrected Hamiltonian H® av transforms as a scalar only if a is independent of E^ since 
E v is the only basic quantity of density weight one. However, the vacuum algebra is first 
class even if a depends on E^. In contrast, when a — 1 (or more generally a/a),a must 
be independent of E v . (The case a = a in vacuum is special because any such correction 
could be absorbed in the lapse function, making the algebra formally first class.) 
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In summary, for corrections a (and a) independent of we have 

{H^JNID^IN*)} = -H^ V [N*N') 

and 

{H£JM),H^[N}} = D grav {a 2 \E*\(En- 2 (MN>-NM>)] (60) 

-G^ia^E^iE^iMN' - NM')rj'\ 

To proceed, we will include matter in the form of dust as it is assumed in LTB models. 
4.2 Dust 

For a full consistency analysis based on the constraint algebra we have to use a dynamical 
formulation of the dust matter source, rather than a phenomenological implementation via 
the dust profile F(x). It is convenient to use a canonical formulation for dust with stress- 
energy tensor T Q/3 = tUJJp as developed in [32]. The dust four- velocity is given by the 
Pfaff form U a = —r, a +WkZ k , a , where as canonical coordinates the dust proper time r and 
comoving dust coordinates Z k with k = 1,2,3 appear. Their respective conjugate momenta 
will be called P and P&. Matter contributions to the diffeomorphism and Hamiltonian 
constraint read 

D dnst [N a ] = J d 3 xN a D a = J d 3 xN a (Pr, a +P k Z k , a ) 
H dnst [N] = J d 3 xN^P^ + q- b D a D b . (61) 

Imposing spherical symmetry and using adapted coordinates $ := Z 1 , Z 2 = Z 3 = (p 
the constraints become 



D dm [N x ] =4tt J dxN x {P t t' + P$&) 



H dust [N] =4n I dxN^P 2 + j§^(PtT' + P^') 2 (62) 

with the remaining canonical pairs 

{r,P T } = {$,P*} = i- 

whose momenta P T and P$ are defined by the relations P = P T sin d (in terms of the P of 
the full 3-dimensional theory) and P$ = —P T W\. 

For non- rotating dust, as must be the case with spherical symmetry, the constraints 
Pfc = can be imposed by requiring that the dust motion be described with respect to the 
frame orthogonal foliation, so that the state does not depend on the frame variables Z k . 
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As a result P$ is usually taken to be zero. However, we will not choose to do so until we 
try to solve the equations of motion. 

From the form of the Hamiltonian (IBTT) in the full theory and q ab = (det E^,)~ 1 E°'E b , 
we can expect quantum corrections P[E X , E v ] from a quantization of inverse triads inside 
the square root: 



Also here, the form of j3 will be restricted by the requirement of anomaly freedom. 

Adding the individual contributions, the diffeomorphism and corrected Hamiltonian 
constraint for the gravity-dust system are D[N X ) = D gr3iV [N x ] + D dust [N x ] and H Q [N] = 
H^ rscv [N] + H® ust [N] . Now, the Poisson bracket for the matter part of the Hamiltonian with 
the diffeomorphism constraint is 



{HZ s tW},D[N x ]} = -HLtW x N>] + J dxNN x 



, dP \E X \ D 2 X 



with D x := P t t' + P^'. The closure of {H Q [N}, D[N X }} consistently imposes the condition 
that a and (3 be independent of E v , upon which 

{H Q [N],D[N X }} = -H Q [N X N'} . 

Finally, 

{H® [N] , H* [M] } = {Hg av [M] , Pg av [N] } + {H^ st [M] , H$ ust [N] } 

= D gIiiV [a 2 \E x \(E ip )- 2 (MN' - NM')] - G griw [a 2 \E x \(E ip )~ 2 (MN' - NM' 
+ D Aust [{3\E x \(E*)- 2 (MN' - NM')} 

gives the relation 

(3 = a 2 . (63) 

Note that the presence of matter makes this consistent deformation of the classical con- 
straint algebra non-trivial: corrections can no longer be absorbed in the lapse function. 
We also point out that a deformation of the constraint algebra is required to implement 
the corrections consistently. This seems to be an interesting difference to a reduced phase 
space quantization which is possible in this class of models based on a deparameterization 
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4.3 LTB-like solutions 



Using the transformation equations (EJ), the quantum corrected Hamiltonian in ADM vari- 
ables © reads 



H Q [N) = ^ y"dxAT 



PlPr , _ , LP? L R 1 2 

a 1 +(2a — a) — a (2a — a) h a 

R v ; 2i? 2 2 v 7 2L 



RR'\' 



+ AnGP T ( 1 + /3— 



./ 2 



and the diffeomorphism constraint is 

D[N X ] = 1 y dxN x {R'P R - LP' L + 47rGP r T / ) 

where we have already used P$ = 0. Paralleling our treatment of LTB-reduced constraints 
our further analysis will be split into two different cases of correction functions. We choose 
to work here in ADM variables instead of triad variables, but of course identical results 
follow using the latter. 



4.3.1 Case a = a 

The equations of motion R = {R,H Q [N] + D[N X ]}, L = {L, H Q [N] + D[N X }}, P R = 
{P R ,H Q [N} + D[N X }}, P L = {P L ,H Q [N} + D[N X ]}, f = {t,H q [N] + D[N X ]} and P T = 
{P T , H Q [N] + D[N X ]} are respectively, 

P L = — (-R + N X R') (64) 



aN 

Pr = ^n (~ Lk ~ Lr + ( nxrl Y) ( 65 ) 

• Ar (PlPr LPl\ Ar da / P L P R LP? L 

P R = -Na - N — H - 

R V R 2 R 3 J dR\ R 2R 2 2 

Ar R'Y Ar , R' / Tl R\ , da RR' 
Na— + N'a N'a— + N' 



L J L V LJ dR L 

-Z*GNP T (l + P-) -JL (66) 

r /2 / /2 



+ AttGN(3P t — \l+(3—) (67) 
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f = N[l+ p^P) 2 + iVV 



Pr 



N(3 



L 2 



-I 2 



+ N X P T 



(68) 
(69) 



To try to find an LTB-like solution, following [16] and the recent [33] , we choose an embed- 
ding by coordinates such that r = t, or equivalently, from (1B"81) . N — 1. Setting iV^ = 0, 
equation ([691 becomes P T = so that P T (x) is a function of the spatial coordinate only. 
Substituting this and fl64l) with fl65|) in the diffeomorphism constraint 



5D_ 

In* 



R!P R - LP' L + AttGP t t' = 



gives the equation 



aL 



or R' /aL = £, with £(x) an arbitrary function of the spatial coordinate only. To make 
contact with the classical LTB solution we may define k by 



(70) 



and obtain the corrected LTB condition 

at 

Note that this is exactly the form of corrections used earlier in (j23|) . which we now have de- 
rived from the corrected constrained system without extra assumptions about the possible 
form of LTB solutions. 

With all this, the Hamiltonian constraint 



a 



PlPr LP* 
2R 2 



R 



L 

2 



R' 2 ( RR! 

+ 



2L 



L 



AnGP T = 



becomes 



RR 2 



a z 



+ R(l-a 2 E 2 ' 



8ttG£P t . 



(71) 



Again, we may define F by the equation 8ttG8,P t = F', so finally the integrated Hamilto- 
nian constraint reads 

RR 2 = a 2 F + Ra 2 (a 2 8, 2 - 1) + c(t) (72) 
and the second order, evolution equation from (|67|) is 



2RR + R 2 = 2(R 2 + a 4 £ 2 )^^ - n~( I - ( r'£- 

d log R 



(73) 
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This is precisely the time derivative of (172"]) . provided c(t) = c = const, and shows consis- 
tency. In the limit a = 1 we recover the classical LTB condition and evolution equation. 
In this limit the integration constant can be absorbed in F so we may set c = here. 

Even though the corrected LTB condition coincides with the one derived earlier, the 
first order and evolution equations are slightly different from the corresponding ones (1391) 
and (H2J) derived before: 



RR 

a 2 



2* ' 



+ R'{1 - a 2 £ 2 ) - 2EE'ahR = F' (74) 



which unlike (1711) is not a spatial derivative, and the second order equation 

2RR + R 2 = 2R 2 ^^- -a 2 (l- a 2 E 2 ) . (75) 
d log R 

While several terms in the two sets of the equations agree, the spatial derivative of (|72p dif- 
fers from flTJD by aterm2££ , a J R(/i-a)+2£ 2 iW, and (J7J) from fl7S} by 2a 4 £ 2 dloga/dlog J R. 
It is interesting to note that the more restrictive new method of this section provides a 
Hamiltonian constraint which can be spatially integrated. 

4.3.2 Case a = 1 

Similarly to the previous case, choosing N — 1 and N x = gives the equations of motion 

P L = -RR (76) 
P R = - (RL)' - (1 - a)LR (77) 

Pl=1 + (2-«)(-^-^] (75 



n 


R'2\ 




2R 2 


2L 2 ) 








da 


)- 




" dR 



LP? P l Pr\ fR'X da (LP 2 L R' 2 



Fr ' R 3 R 2 ) { L ) + dR { 2R 2 H ^ 



2 2L 



+ (l-«)^-((l-«)f)'. (79) 

Substituting the first two equations in the diffeomorphism constraint gives 

R\' , .RR 
1 — a) 



n L I v ' RL 
Inserting the ansatz ( flTj) 

£/ 

gives the same equation ( TT81) for the function f(E x ) as in the previous treatment: 

= /(l - «) • (80) 
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Substituting equations (!76|) and (ITT)) into the Hamiltonian constraint gives 

(2 - a)RR 2 + 2RRR - 2RR'Rj + aR - aRE 2 f 2 - R{1 2 f 2 )' = 8irGEfP T 

or using (l80i) and defining F' = 8ttG£P t 

aR'R 2 + 2RRR' + aR(l - 8, 2 f) - R(E. 2 f 2 )' = fF' (81) 

which reproduces the classical equation in the limit a = f = 1. 
The second order equation from (ITS]) is 

2RR + aR 2 = -a(l + £ 2 / 2 ) + 2£ 2 / 2 . (82) 

Also here, compared to ( 1431) and ( )44l) as obtained earlier: 

a^R'i? 2 + 2RRR + aR(l - £ 2 / 2 ) - R{£ 2 )' ' fl = fF' 

2RR + aR 2 = -a(l - £ 2 / 2 ) , 

the structure of the resulting equations remains similar up to a few extra terms. 

As in sections 13.51 and I3T61 we may interpret the effects of correction terms using effective 
densities and energy functions. 

5 Possibility of singularity resolution through bounces 

We can now use the different sets of consistent equations to analyze properties of gravita- 
tional collapse, such as the formation of singularities. Classically we have the first order 
equation 

R 2 = k + - . (83) 
R 

Compared with the marginal case where k = 0, there exists the possibility that R = even 
for positive mass functions F. However, to conclude whether there is a bounce or not we 
also need to look at the evolution equation 

2RR + R 2 = k (84) 

and see whether we can have R > in addition to R = 0. From the first equation we see 
that R = implies k + F/R = and with both F and R positive we get k < 0; bounces 
would be possible only for negative k. On the other hand, for R = the second order 
equation implies 2RR = k. Since R > and k < we conclude that R < and thus there 
is no bounce classically. 

We would like to proceed in a similar manner for the quantum corrected case also. 
However, there the first order equation in time contains terms with spatial derivative as 
well (which cannot be integrated out in all cases). Therefore the analysis for the possibility 
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of a bounce cannot necessarily be done as easily as for the classical case, a feature clearly 
related to the fact that we are dealing with inhomogeneous models. Furthermore, because 
of the inhomogeneous nature of the problem a bounce also makes the analysis difficult by 
the fact that after the bounce there will be the possibility of shell crossing (unless shells 
with larger values of x bounce at larger values of R). 

We note that fl39|) and (|4"3l) imply, as for the classical case, that R = is possible in 
both versions of inverse triad corrections. Whether this corresponds to a bounce is what 
we want to analyze. We start by putting R = in fl39l) which gives 



2 2' 

a — na 



R' - n'ahR = F' . (85) 



Similarly, the evolution equation (j4"2l becomes 

2RR = -a 2 (l -a 2 - na 2 ) . (86) 
Using (j55|) on the right of the above equation we have 

« ^F' + n'ahR . . 

H = - a ^Ri^- (87) 

For a bounce, R > which implies that — a 2 (F' + hz'ahR)/2RR' should be greater than 
zero. We need to check whether this condition can be satisfied in the non-marginal case 
where, as mentioned before, there are two possibilities k > and — 1 < k < 0. In what 
follows we will assume that R' > 0, which locally around a potential bounce is a valid 
assumption even though a collapsing shell with x = Xi, say, would start expanding after it 
experiences a bounce: when R(t,Xi) = the radius of this particular shell is not changing 
with time. The assumption then essentially implies that a shell with x — x^ > xi, which 
may still be contracting, does not immediately catch up with the x\ shell. With this 
assumption and because R and a are positive, the condition that (IH7I) be positive becomes 

k' < -F'/ahR (88) 

as a condition on k' which needs to be satisfied for a bounce, implying in particular that 
for a bounce k' has to be negative. Whether fl88l) can be satisfied generically is not clear 
and must be determined from a numerical analysis of the equations. 

For the first version of the inverse triad correction, R = in (T4"5|) gives 

a(l-f 2 -Kf 2 )R'-K f flR = fF f . (89) 

The evolution equation ( 1441) becomes 

2RR = -a(l - f - Kf 2 ) (90) 

and using (1591 on the right we get 

F' + k'IR , . 

R=-f — 1 . 91 

J 2RR> v ; 
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This should be greater than zero for a bounce. Again because of the presence of spatial 
derivatives in the above expressions it is difficult to say whether there is a bounce in general 
and whether or not the singularity can be avoided. 

With (|72l) we have a corrected equation which can be spatially integrated, allowing an 
analysis similar to the classical one. The condition R = at a bounce implies 



1 — a 2 £ 2 1 — a 2 — K,a 2 
which is positive provided £ < 1/a. The second derivative 

2RRU = 2a 4 4^£ 2 - « 2 (1 - a 2 £ 2 ) 
lR -° dlogi? v ; 

can be positive under this condition only if the derivative dloga/dlogi? is sufficiently 
positive. This is not the case in semiclassical regimes (for geometries to the right of the 
peak in Fig. [I]), where bounces are thus prohibited. The correction function is increasing to 
the left of the peak, which is a regime with strong quantum geometry corrections. Since we 
have not included all quantum corrections, a conclusion of a bounce in this regime would 
be unreliable. The only semiclassical option for a bounce is to have a geometry above 
the peak of inverse triad corrections, but have decreasing patch sizes A which appear as 
the argument of a. Thus, the number N of patches would have to increase sufficiently 
rapidly. In this regime, we have a > 1 and thus £ < 1 by our condition for R > 0. Such 
a bounce would thus be possible only for k < 0. (As the argument shows, without lattice 
refinement a bounce from inverse triad corrections would at best be possible only in the 
strong quantum regime.) 

While bounces seem possible in the present situation, they cannot be considered generic. 
They require a regime where the patch number is increasing sufficiently rapidly in such a 
way that the patch size decreases. Since in the discrete quantum geometry of loop quantum 
gravity the patch size has a positive lower bound, the patch number of an orbit of fixed 
size cannot increase arbitrarily. Tuning then seems required to have the right behavior just 
when a shell is about to bounce. 



5.1 Near center analysis regarding bounces 

On the basis of a general analysis it seems difficult to conclude whether quantum corrections 
generically resolve the singularity in non-marginal LTB models through bounces. The 
simplest possibilty seems to be one where the central shell is prevented from becoming 
singular because of a bounce. For almost complete collapse, we should expect the relevant 
regime to be one of small R. In this case the subsequent study of the outer shells will 
become difficult due to possible shell crossings, but presumably these outer shells will not 
become singular either. We therefore now proceed to a near center analysis. 

As in the marginal case of [3], we use techniques similar to those in [M] and assume 
that near the center of the dust cloud we can expand R(t, x) as 

R(t,x) = R 1 (t)x + R 2 (t)x 2 + --- (92) 
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For the classical collapse the mass function can be expanded as 

F(x) = F 3 x 3 + F 4 x 4 + • • • (93) 

Substituting the expansion for R(t, x) and for F(x) in the classical first order equation 
R 2 R = kR + F, we find that the lowest order term on the left (as also the second term on 
the right) of the above equation goes as x 3 . This suggests that the series expansion for the 
energy function k(x) should be 

k(x) = k 2 x 2 + n 3 x 3 H (94) 

It also implies that at x — 0, the center of the cloud, k(x) = and therefore for the case 
where k(x) > outside x = 0, k 2 should be greater than zero. On the other hand, for 
— 1 < k(x) < 0, k 2 should be less than zero. However, if we consider our effective k as in 
(j54p then we can have lower order terms in k. Since the lowest order term in a is of order 
x 3 and that in / is of order x we can have the lowest order term in k behave as x~ 3 . This 
would then imply that at the center k blows up whereas classically for negative k we have 
the condition — 1 < k < 0. 

With this caveat we now consider ( 1891) to lowest order after using various series expan- 
sions: 2 

Cl R 3 x 3 (l - c\R\x 2 - C2K ~^ Rl ^j R i + 3c 2 c 3 K_ 3 Rtx 2 = 2c 2 F 2 R l x 2 . (95) 

Here c\, c 2 , c 3 and k,- 3 are the coefficients of the lowest order terms in the expansion of a, 
f, I and k, respectively. Here we assume the orbital vertex number N to be nearly constant 
around the center. (This gives rise to the strongest effect from inverse triad corrections 
and allows direct comparisons with the matching results from [13J.) The lowest order term 
in F is F 2 instead of F 3 because of the effective density correction. For x ~ the first two 
terms in the parenthesis on the left can be ignored compared to the third term and thus 
the above equation implies that R\ = for 

_ 2c 2 F 2 

K - 3 (3c 2 c 3 - Cl c 2 )Rl ■ l9bj 

It turns out that 3C2C3 — C\c\ < implying that R\ can be zero only for k_3 < as in 
the classical case. If we now look at (1901) near the center we find that the condition for a 
bounce is 

p 

K -> > 3^ < 97 > 

which means that, as in the classical should be positive implying that we do not 

have a bounce for the central shell. 

For version two the expression for the effective density (1531) and the effective energy 
function (j55p . respectively, imply that to lowest order the mass function can behave as x~ 3 
and the energy function as x~ 10 . To look at the possibility of a bounce at the center of the 
cloud consider fl85l) to lowest order: 

D ^2p7f! g-lo4gl 1 10ft-10CiC 4 itf 3F_ 3 

Ri-<h R i x -4 — + -4 = — ^r- ( 98 ) 
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Here C4, _F_ 3 and k_io are the coefficients in the expansion of h, F and k respectively. 
Ignoring the first two terms for x ~ 0, the above equation implies 

3F 

K ~ 10 = _ 7TFi ~ 3 2n P 7 (") 

(10cic 4 — cf)R{ 

as the condition for R\ = 0. We note that the denominator here is positive implying that 
K-io > if F_ 3 < (negative mass function near the center) and k_io < if F_ 3 > 
(positive but decreasing mass function near the center). None of these behaviors could 
occur classically; either negative total energy (F < 0) or a negative density (F' < 0) would 
be required. To see if the above condition implies a bounce we use (1861) to lowest order 
and find 

3F 

K -» > -lo^k (100) 

as the condition for getting a bounce. This means that for F_ 3 < 0, k_ 10 has to be positive 
(in agreement with the condition for Ri = found above) implying that a bounce for 
the central shell is possible if the above inequality is satisfied. For F_ 3 > 0, k_ 10 has to 
be greater than a negative number and thus if it is positive then a bounce again seems 
possible. 



6 Collapse behavior near the center 

Using (jS2D, dnSD and (1541) we see that to lowest order in x (which is x 3 ) the classical equation 
R 2 R = kR + F implies 

dR 1 



dt 



K 2 + 



F3 
Rl 



This has the solution (choosing the minus sign which corresponds to collapse) 



(101) 



iW«2 + § F 3 



K2 



+ — ^log(F 3 + 2k 2 R, + 2^ x /k 2 + - f ^R 1 



2k 



Ri 



(102) 



for k 2 > and 



«2 



«2 



^3 

+ 771 — arctan 

2 «2 7 



f -| K2 |(2| K2 |i2 1 -F 3 
2^2T(«: 2j Ri-F 3 ) 



(103) 



for k 2 < 0. 

We now proceed with a similar analysis for the first version of the inverse triad correc- 
tion. Near the center the various quantities (a, f, I) behave as 



a 



3/2 



R\x 3 , f 



,l-7r/2 



ill 



-R\x , I 



1 / 2 \ 3/2 /Se 1 --/ 2 



i^x 4 (104) 
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where the way the near center behavior for / has been determined is described around 
Eq. (1321) . In what follows we will denote the constants appearing in the expansion of 
(a, /, I) by (ci, c 2 , c 3 ) respectively Thus substituting the series expansions in ()43l) we get 

ci(l - c\R\x 2 - k 2 c 2 2 R\x a )R{x 3 + c x R\R\x h + 2R 2 1 R 1 x 2 - 2k 2 c 2 c 3 R\x 1 = 3c 2 F 3 R 1 x 3 (105) 

We now consider three different possibilities. 

6.1 Case I: No correction to the expansion of F and k 

If we work with (11051) directly then we find that the lowest order term on the left hand side 
is 2R\Rix 2 and the lowest order term on the right hand side goes as x 3 implying Ri = 0. 

6.2 Case II: Modification to the expansion of F and no modifi- 
cation to K 

However because of the presence of an extra factor of x on the RHS we can start the 
expansion of F with the leading term behaving as x 2 . In this case the lowest order term 
on both the LHS and RHS are of order x 2 and we get 

2R\R l x 2 = 2c 2 F 2 R l x 2 (106) 

which has the solution (for collapsing dust cloud) 

R 1 (t) = l-^~F 2 (t-t ) (107) 

where we choose the initial condition R\{t = t ) = 1. We see that the central singularity 
Ri = forms in a finite time t s — (1 + \J c 2 F 2 t^) j ^/ c 2 F 2 . 



6.3 Case III: Modifications to the series expansion of F and k, 

There is a third option which is suggested by the possibility of a corrected energy function as 
discussed earlier. If we consider this correction then the lowest order term in the expansion 
of ft goes as K-3/x 3 . With this included the matching of lowest order terms in (11051) gives 



- k^ Ci c 2 2 R\ + 2R\ + Zk^c 2 CzR\ = 2c 2 F 2 (108) 



and thus 



with the solution 



dt = ^ (109) 

\/2c 2 F 2 + (K_ 3 d4 - 3/c_ 3 c 2 c 3 )i2f 



[K r> 1 1 1 (K-3Cic|-3K-3C2C3)-Rf 77/1 1 6 (k_3Cic2-3k_3C2C 3 )_R5 x 
V2^1-\/l + ^ 2M5>2>5> 2^F 2 ) , , 1ln x 

t = , h c (HO) 

a/2c 2 F 2 + (k_ 3 cic| - 3k^ 3 c 2 c 3 )RI 

where 2 Fi(a, b; c; x) is the hypergeometric function and where Co is constant of integration. 
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7 Homogeneous limit and Minkowski space 



In the classical case, the first order equation is R 2 R = k(x)R + F(x) with corresponding 
expression for the mass function F' = 8iiGR 2 R'e. This allows isotropic space-times as 
special solutions. We use the ansatz R(t, x) = a(t)x with the condition that at time t = 0, 
a(0) = a and assume that the density e = eo is a constant and that the energy function 
goes as k = —kx 2 where A; is a constant. When used in the first order equation, this gives 

.2 . . 87rGe , . 

a a = —ka H (111; 

3 

which is the Friedmann equation. 

We would now like to see if we can get a Friedmann-like solution within the LTB 
class with inverse triad corrections included. This would indicate whether there can be an 
effective geometry of the classical homogeneous form, although the notion of homogeneity 
itself might change on a quantum space-time. Using the ansatz for R and the assumed 
form for the mass function and the energy function we find that for the first version of the 
inverse triad correction (|43[) we get 

ail — f 2 ) , „ 9 n aa 2 ax 2 2aa 2 x 2 n . , 9 „ _ 9 + 
v - 1 a + kafax 2 + — - — + — - — + 2klax 2 = 8nGe x . (112) 

Although the resulting expression is not as simple as in the marginal case, we can see that 
Friedmann like solutions are not possible since x 2 does not cancel from the first term while 
a is allowed to depend only on t. 

The second version of the inverse triad correction (|39j) gives 

(1 - a 2 + kc?x 2 )a + 3 ° - 2a f f a + 2kahax 2 = 8ixGe Q x 2 . (113) 

a 2 cr a(ax) 

Again the Friedmann solution is prohibited. Since the first term, which spoils the homo- 
geneous limit, is the same for an analysis based on (172"]) . there is no homogeneous limit 
in that case, either. One can also see from the first term that no other x-dependent k, 
which might implement quantum corrections to the notion of homogeneity, can resolve the 
non-existence of homogeneous effective geometries subject to our equations. 

The homogeneous limit, as a special case, would also include Minkowski space as the 
vacuum solution. Classically the first order equation for F = and k = 0, implies that 
R = R(x) and we recover the Minkowski spacetime. However from (H3l we see that after 
choosing the mass function and the energy function equal to zero the equation becomes 

ail - f)R' + aR 2 R' + 2RR'R = (114) 

which, due to the presence of the first term, implies that R will be dependent on both 
(t,x). Even though the equation has the correct classical limit, it is not straight-forward 
to see how the time dependence of R should disappear in the Minkowski limit. 
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For the second version (139]) gives 



... .1, „. R 2 R' 2RRR' 2RR 2 R' da 

1 - a 2 )R' + + — — = 115 

or a z cr dR 

which again implies that even though in the classical limit we do have a Minkowski solution, 
there is still time dependence in R in corrected solutions. 

Strong corrections are suggested at small values of the argument of a, which, given 
that R determines that value, may seem unacceptable because the center in Minkowski 
space is not physically distinguished. However, the radius R and thus the center is directly 
relevant for the size of corrections only if there is no lattice refinement in which case the only 
parameter which a depends on is R. The primary argument of a is, however, the size A 
of discrete patches rather than R 2 of whole spherical orbits. With a non-trivial refinement 
scheme, a(i? 2 /}sf) will also depend on the number of vertices per orbit, which for a good 
semiclassical state must provide a more uniform distribution of quantum corrections not 
distinguishing a center: N must be small when R 2 is small. If discrete patch sizes on all 
orbits are nearly similar, quantum corrections are uniform and do not distinguish a center. 
A detailed discussion would go beyond the scope of this paper, but one can already see the 
crucial role played by lattice refinement for the correct semiclassical limit. 



8 Matching of the loop corrected LTB interior with 
generalized Vaidya exterior 

To verify the corrected mass formulas from a different perspective, it is instructive to 
match our corrected models of dust collapse to radiative generalized Vaidya solutions. The 
interior metric for an inverse triad corrected non- marginal LTB model can be written in 
the form 

ds 2 = -dt 2 + - W \ dx 2 + R 2 dQ 2 (116) 
(1 + n)q 2 

where q takes the appropriate form depending on whether the correction function corre- 
sponds to first version or the second version. Following [13J, we match this interior solution 
with the generalized Vaidya metric 

ds 2 = - (l - dv 2 + 2dvd X + X 2 dn 2 (117) 

where M is allowed to be a function of v as well as \. 

Let the boundary of the dust cloud be at x = x\>- Here, the exterior coordinates will 
then be functions v = V\,{t), x = Xb(t) of the interior coordinate t such that 

dv = t>bdt , d% = Xbdt . (118) 

With this, the metric induced from the exterior on the matching surface becomes 

ds 2 = - ((l - ™] v 2 h - 2v hX ^ dt 2 + xldn 2 (119) 
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Matching the first fundamental form at the boundary (where dx = in the interior coor- 
dinates) we get 

- dt 2 + Rldtf = -((1-^%- 2v h x^ dt 2 + xldtf (120) 

where i?b is the value of the area radius at the boundary. Equating the coefficients of dt 2 
and dQ 2 , respectively, on the two sides implies 

Xb = Rb 

2M\ „ , , 

)vl-2v hX b = 1. (121) 

Xb / 

We now match the second fundamental forms (extrinsic curvature) at the boundary. For 
this we need the normal to the boundary in the interior as well as the exterior coordinates. 
In the interior we have dx = which implies that the normal = (0, a, 0, 0) where a 
is a constant fixed by the normalization ^fn* nj, = 1 (the label i standing for 'interior'). 
This gives a = R' h /%\/l + k and the normal in the interior is = (0, R^/q^l + n, 0, 0). 
Evaluating the 99 component of extrinsic curvature K l v = nj^ we find 



JC m = q h R h VT+K (122) 
From (11181) . at the boundary, we have dv — v^,dt = 0, dt = ^ and thus 

Xb 

- Xbdv + v h d X = (123) 

In terms of the exterior coordinates the normal to the boundary is given by = (— cxb, cv^, 0, 
where c is a constant fixed again by the requirement g^n e ^n e v = 1. Using the inverse of the 
metric in the exterior we find that c = 1. The 99 component of the extrinsic curvature in 
the exterior is 

/ 2M\ 

K e ee = -XbXb + Xb 1 }Vb. (124) 



Xb 

Equating f)122p and (I124p and using (11211) to simplify the resulting expression we obtain 



/ • ( 2M\ 

q h Vl + n = -R h + M - — 



Vb (125) 



implying 



q h y/T+K + R h 

Vh = f3lM ' ^ 126 ^ 

Using this and its square, in the second line of (I12ip we have 

2M = (I - ql(l + k) + Rl)R h . (127) 
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At this stage, we can note that the expression for the Vadiya mass resulting from matching 
is the same as the one obtained for the effective Misner-Sharp mass in (|50|) and (152|) . 

The exterior region can contain trapped surfaces when the condition 2M = \ is satisfied. 
Using the first expression in (1121R along with (11 2 7ft this implies that the boundary will be 
trapped when (1 — q^(l + k) + RfyRb = Rb, or 



where the negative sign has been chosen since we have a collapsing scenario. This formula 
shows that the horizon condition is corrected by % simply by multiplying 1 + k, 2 as it is 
suggested by the effective metric (11161) . 



We have extended the treatment of [3] to non-marginal models, where additional corrections 
from spin connection terms arise. With these additional terms the original derivation 
appears more arbitrary, which led us to provide an independent derivation of equations 
corrected by the treatment of inverse triads in loop quantum gravity. In this new derivation, 
anomaly-freedom is implemented first and LTB conditions are imposed afterwards to select 
a special class of solutions. The structure of the resulting equations is very similar in 
both derivations, showing the robustness. By the alternative method, which is much less 
arbitrary than the one extended from [3J to non-marginal models, we thus show that the 
more phenomenological treatment of corrections used in [3] is reliable. In details, however, 
the resulting equations do differ which is always possible due to quantization ambiguities. 
The effects analyzed in this paper do not appear to depend sensitively on the method, 
but further analysis may well provide restrictions on acceptable equations, and thus on 
quantization ambiguities. 

Our analysis in this paper has been done for inverse triad corrections, while holonomy 
corrections, which to some degree were treated in the marginal case of [3], are technically 
more involved. Already for inverse triad corrections, the extension provided here is an 
interesting step in the analysis of inhomogeneous collapse and singularities. Comparing 
with homogeneous models and matching results of [13] would suggest easy resolutions 
of singularities by bounces. Marginal models were not entirely conclusive in this regard 
since their homogeneous analog is that of a spatially flat Friedmann-Robertson-Walker 
model which under inverse triad corrections gives rise to bounces only with a negative 
matter potential |35j. In the collapse analysis, however, we have used positivity conditions 
for the mass function which indicate that bounces in marginal LTB models with inverse 
triad corrections should not be expected. For non-marginal models, on the other hand, 
homogeneous special solutions with positive spatial curvature exist, which do show bounces 
with inverse triad corrections and positive matter terms [12] . One would thus expect non- 
marginal models to result in bounces much more easily than marginal ones do. 

This, however, is not the case: we mostly confirm the results found in marginal mod- 
els where (i) bounces are not obvious and (ii) a homogeneous limit of quantum corrected 
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solutions may not even exist. As for the first property, bounces seem somewhat easier to 
achieve than in marginal models, but in contrast to the expectation turn out to be hard to 
realize generically. Moreover, a complete analysis would have to involve an investigation of 
shell-crossing singularities which can be involved even classically. (See [3] for more discus- 
sions on this in marginal corrected models.) As for property (ii) about the homogeneous 
limit, one can evade ruling out a homogeneous limit at the dynamical level only if one 
assumes a quantum notion of symmetry which would imply effective isotropic space-time 
metrics different from classical Friedmann-Robertson-Walker models. This may well be 
expected, as indeed the deformed constraint algebra ( |60j) shows that there is a corrected 
quantum space-time structure. It would be interesting to see how this influences space-time 
symmetries. 

Finally, several issues discussed here involved the role of lattice refinement for the semi- 
classical limit. As treatable models between homogeneous ones and the full theory, LTB 
models turn out to be quite instructive. This should also be expected for an implementa- 
tion at the state (rather than phenomenological) level which was started in [3] for marginal 
models but which we have not attempted here for non-marginal models. 
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